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We propose a method to obtain adiabatic internuclear potentials via energy variation 
with the intercluster-distance constraint. The adiabatic 16 0+ 16 ' 18 potentials obtained 
, by the proposed method are applied to investigate the effects of valence neutrons in 

16 0+ 18 sub-barrier fusions. Sub-barrier fusion cross sections of 16 0+ 18 are enhanced 
more compared to those of 16 0+ 16 because of alignment of valence neutrons in 18 0. 
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■ 1. Introduction 

Recent theoretical and experimental studies have revealed various exotic structures such 

■ as neutron halo and skin structures, where valence neutrons play important roles. Nuclear 
reactions at low incident energy are efficient tools to investigate the effects of valence neutrons 
in nuclear dynamics. For instance, during tunneling through the Coulomb barrier, in sub- 
barrier nuclear fusion reactions, valence neutrons probably affect excitation of nuclei by 
influencing the alignment, polarization, and vibration. To investigate low-energy nuclear 
reactions theoretically, internuclear potentials that work well in low-incident energy are 
required. 

There are two types of internuclear potentials for nuclear reactions; adiabatic and sudden 
potentials. To obtain adiabatic potentials in mean-field approaches, two methods have been 
proposed: applications of the time-dependent Hartree-Fock (HF) method with a density 
constraint [14] and the HF method with an internuclear-distance constraint in a symmetric 
form[17, 18]. In the former method, energy is minimized using the constraint on the density 
distributions obtained by the time-dependent HF calculations. Since the density is deter- 
mined by time-dependent HF at incident energies greater than the Coulomb barrier, the 
distortion effects of colliding nuclei associated with density changes at low incident energy 
(i.e., below the Coulomb barrier) can be insufficient. In the latter method, energy variation 
with a symmetric constraint on the internuclear distances has been applied to symmetric 
systems but not to asymmetric systems. Double-folding potential, which is a type of sud- 
den potential, has also been used to study nuclear reactions; fusion cross sections near the 
Coulomb barrier are described by channel-coupling (CC) calculations with double-folding 
potentials [4]. Recently, a repulsive core potential has been suggested to account for deep 
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sub-barrier fusion cross sections [8]. However, the fundamental origin of the phenomenological 
repulsive effect has not been clarified yet. 

In sub-barrier nuclear fusion, adiabatic potentials are expected to work, which contain 
excitation effects of colliding nuclei. Adiabatic potentials treating excitation effects of col- 
liding nuclei contains CC effects approximately [2, 3, 9, 15]. To investigate the excitation 
effects of valence neutrons on fusion cross sections, comparison of 16 0+ 16 ' 18 sub-barrier 
fusions are feasible because 18 possesses two neutrons more than 16 0. The sub-barrier fusion 
cross sections of 16 0+ 18 are enhanced more compared with those of 16 0+ 16 although 
the enhancement is smaller than that in heavier systems such as Ca isotopes[13]. The CC 
model shows that most of the enhancement is attributable to the excitation of the low 2 + 
state in 18 0[16]. In potential model interpretation, it means that the excess neutrons result 
in an effectively thinner or lower Coulomb barrier for 16 0+ 18 than that for 16 0+ 16 0. 
The CC model is a popular to analyze fusion reactions around the Coulomb barrier and 
gives reasonable results [2]. However, most of internuclear potentials used in the CC model 
are phenomenological ones, and they are not based on microscopic frameworks where the 
antisymmetrization effects between colliding nuclei are taken into account. 

This paper aims to propose a method to obtain adiabatic internuclear potentials in a 
full-microscopic framework. In the process, the energy is minimized while constraining the 
internuclear distances via the deformed-basis antisymmetrized molecular dynamics (AMD) 
model, which can be easily applied to both asymmetric and symmetric systems. The adi- 
abatic 16 0+ 18 potentials obtained by the proposed method show that the excitation of 
valence neutrons in 18 O reduces the thickness of the Coulomb barriers and enhances the sub- 
barrier fusion cross sections of 16 0+ 18 O. The excitation effects are analyzed via comparison 
of cross sections obtained by adiabatic and sudden potentials. In general, the enhancement 
effects of sub-barrier fusion cross sections are larger in heavier systems but we choose the 
light systems 16 0+ 16 ' 18 0, due to the numerical cost for obtaining the internuclear potentials. 

In Sec. 2, we explain the framework to obtain internuclear potentials and fusion cross 
sections. In Sec. 3, we present the internuclear potentials and fusion cross sections. In Sec. 4, 
we discuss the role of valence neutrons in 18 to enhance sub-barrier fusion cross sections. 
Finally, conclusions are given in Sec. 5. 

2. Framework 

We propose a method to obtain adiabatic potentials via energy variation with the 
intercluster-distance constraint in the AMD framework[12]. In the present study, we use 
the deformed-basis AMD framework[7]. A form of the deformed-basis AMD wave function 
|<&), Slater determinant of Gaussian wave packets, is described as 



parts, respectively. K is a real 3x3 matrix that denotes the width of Gaussian wave packets, 
which is common to all nucleons, and Zj is a complex vector that denotes a centroid of a 
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Gaussian wave packet in phase space. The wave function is set as the expectation values of 
the position of center of mass and the total momentum are zero. A wave function |<&Ci-C 2 )> 
having a dinuclear structure comprising nuclei Ci and C2, is defined as 

|$ Cl -C 2 )=i(|$ Cl )®|$C 2 )), (4) 

where |<3?c t ) is a direct product, which is not antisymmetrized, of single-particle wave 
functions with proton and neutron numbers corresponding to the nucleus C«. 

The internuclear distance R is defined by the density distribution of a total system. Suppose 
the centers of mass of nuclei Ci and C2 are located on the z-axis with z < and z > 0, 
respectively. Boundary planes of nuclei Ci and C2 for protons and neutrons are denoted by 
z = z p and z = z n , respectively, and are defined as 



dz II dxdy p p (r) = Z 1: (5) 



dz J J dxdy p n (r) = N\. (6) 

Here p p (r) and p n (r) denote proton and neutron densities, and Z\ and N\ denote proton 
and neutron numbers of nucleus Ci, respectively. The internuclear distance R is defined by 
the positions Ri and R2 of the centers of mass of nuclei Ci and C2, respectively, as 

R = |R.2 — R-i I , (7) 
R. . Z ' R "^ R '"', (8) 

a?'"' = 11 i*&y !'" " iz 1 p *-" {t) < <9) 

J J J —00 

R?' n) = J J dxdy J^dz r p Ptn (r), (10) 

where Zj, iVj, and Ai denote the proton, neutron, and mass number of a nucleus Q, 
respectively. 

To obtain the adiabatic potentials V^d, we optimize the dinuclear wave function while 
constraining the C1-C2 distance using the d-constraint AMD method[12]. That is, the energy 
is minimized while constraining the distance parameter d between the centers of mass of the 
wave packets of nuclei Ci and C2 according to 

5 [(^% 2 ;d\H'\^% 2 ;d) + V cnst (d)] = 0, (11) 
H> = f + V N + V c -2T G , (12) 

where T, Vn, Vq and Tq are the kinetic energy, effective nuclear interaction, Coulomb inter- 
action, and kinetic energy of the center-of-mass motion of the total system, respectively. 
V cns t(d) denotes a parabolic constraint potential for the internuclear distance d defined by 
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using the set of single-particle wave functions |3>Ci) : 

V cnst {d) = v cnst \ (^cJRoiil^) - ($c 2 |RG 2 i|$c 2 > 



(13) 
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where v cns t denotes a sufficiently large number. Details of the constraint potential are 
reported in Ref.[12]. In d > 6 fm, the distance R defined by the density distribution agrees 
with d defined by the centers of mass of subsystems for the 16 0+ 16 system[ll]. Adiabatic 
potentials reflect structural changes with respect to internuclear distances. By using the 
optimized wave function |$q P ^ c an adiabatic potential V a d(R) is defined as 

V ad (R(d)) = (^% 2 ;d\H'\^% 2 ;d) - (E Clgs + E C2gs ), (15) 

where Ec lgs + Ec 2gs denotes a summation of ground-state energies of nuclei Ci and C2 
obtained by varying the energy for isolated systems Ci and C2 in case of common width 
matrices K for the wave functions of Ci and C2. Wave functions obtained by the energy 
variation for a summation of energies of nuclei Ci and C2 are denoted as |$Jf s) ) (i = 1,2). 

To analyze excitation effects of colliding nuclei, we also define the sudden potentials 
V su d(R). We use dinuclear wave functions |<&(? S 1 C2 ; R) defined by ground-state wave functions 
of nuclei Ci and C2. The wave functions \&q^_ c ;R) are defined by shifting the ground-state 
wave functions |$[f s ' ) ) (i = 1,2) to a certain position such that the internuclear distance is 
equal to R, and the total wave function is antisymmetrized. Thus, the structures of nuclei 
Ci and C2 are frozen, except for the effects of antisymmetrization between nuclei Ci and 
C2. Next, we define the sudden potential as, 

V sud (R) = ($ { g?_ C2 ;R\H'\^_ C2 -,R) - (E Clgs + E C2gs ). (16) 

For the 16 0+ 18 system, since the ground-state wave function of 18 is deformed, the 
orientation of 18 is averaged to obtain the sudden potential: 

v sud (R) = -^J v^_ C2 (R,n)dn, (17) 



v^_ C2 (R,n) = (^% 2 ;RMH'\&£% 2 -,R,n) - (E Clgs + E C2gs ). (18) 

For practical purposes, integration is achieved by averaging the direction of the shift in 
the position of nuclei Ci and C2. 

The data points for the potentials are calculated at intervals of approximately 0.5 fm and 
are interpolated by spline curves to obtain potentials as functions of R. 

To obtain the present internuclear potentials, we adopt a distinct treatment of subtracting 
2Tq instead of Tq from the Hamiltonian [Eq. (12)]. The 2Tq is subtracted to eliminate the 
kinetic energy of the internuclear motion T re i and that of the center-of-mass motion Tq of 
the total system. The expectation value Tj of the kinetic energy of the center-of-mass motion 
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of nucleus Cj is separated into the classical part = 2 Am ano - the other part Tqi as 

Ti = Tf ) + T GU (19) 

here Pi is the expectation value of the total momentum of the nucleus Cj. When two nuclei 
Ci and C2 are well separated, Tqi and Tq 2 are functions of K as 

r G1 = T 02 = ^trQ«KK), (20) 
which equals to Tq. For the adiabatic condition Pi = P2 = 0, wave functions of subsystems 

(cl) 

are set as =0 in obtaining sudden potentials. In obtaining adiabatic potentials, A 
term is used as a constraint potential for the adiabatic condition in the 

(cl) 

energy variation, and resultant value of T> is small and negligible. It gives 

Ti = T 2 = Tq (21) 

in both cases of adiabatic and sudden potentials. Using the equations, the 2Tq term is 
written as a summation of the Tq and T TC \ as 

2T G = T 1 +T 2 = T g + T rel , (22) 

and the expectation value of H' is written as 

(H>) = (f) + (V N ) + (V c ) -Tq- T rel . (23) 

Although the Tqi and Tq2 values deviate from Tq in the overlap region because of anti- 
symmetrization, this effect is small in the barrier region because overlap between nuclei Ci 
and C2 is small in the region. Hence, 2Tq is subtracted in the present calculations. By these 
definitions, the internuclear potentials indicate Coulomb potentials in large R region. 

The Modified Volkov No.l case 1[1] and Gogny D1S (D1S) interactions are used as effective 
nuclear interactions Vn- In the Modified Volkov No.l interaction, a three-body contact term 
is replaced with a density-dependent two-body term, and a spin-orbit term of the D1S 
interaction is added to adjust the threshold energy of 32 S to that of 16 0+ 16 (MV1'). 

To obtain fusion cross sections, we use the potential model code POTFUS3 provided by 
Hagino et ai[5]. The POTFUS3 directly integrates second order differential equations using 
the modified Numerov method to solve the Schrodinger equation. Inside the Coulomb barrier, 
the incoming wave boundary conditions that there are only incoming waves at v — ^min is 
employed. The r m j n is set to 6 fm in the present calculations. 

3. Results 

Figure 1 shows the adiabatic and sudden 16 0+ 16 ' 18 potentials calculated for the MV1' and 
D1S interactions. The (T^ + T 2 {cl) ) value is less than 5 x 10" 2 MeV for R > 4.5 fm region. 
The two interactions result in qualitatively similar internuclear potentials. Each sudden 
potential has a structural repulsive core[10] in the R < 5 fm region because of the Pauli 
blocking. In the R > 6 f m region, the adiabatic and sudden 16 0+ 16 potentials are similar 
to each other. Both potentials show barrier tops at almost the same internuclear distances, 
and the shape of the potential curves is also similar within the barrier top. The sudden 
16 0+ 18 potential is similar to the adiabatic and sudden 16 0+ 16 potentials. However, 
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Fig. 1 Adiabatic and sudden 6 0+ ' O potentials calculated as functions of the inter- 
nuclear distance for (a) the MV1' and (b) D1S interactions. Solid and dotted lines represent 
sudden and adiabatic 16 0+ 16 potentials, respectively; short-dashed and dot-dashed lines 
represent sudden and adiabatic 16 0+ 18 potentials, respectively; and the long-dashed line 
represents the point Coulomb potential. 




the adiabatic O+ O potential is lower than other potentials. Due to the lower nuclear 
potential, the barrier tops of the adiabatic 16 0+ 18 potentials occur at larger internuclear 
distances. The difference between the sudden and adiabatic potentials indicates that the 
effect of excitation is large in the 16 0+ 18 system, whereas it is small in the 16 0+ 16 
system. In both potentials, the MVl' interaction gives lower potentials as compared with 
those given by the D1S interaction. 

Figure 2 shows the density distribution of the proton and neutron parts in the 16 0+ 18 
wave functions at R = 7.5 fm. This distribution is obtained by varying the energy with a 
constraint on the internuclear distance for the MVl' interaction. The density distribution 
of the two valence neutrons, which is defined by subtracting the proton density from the 
neutron density assuming density distributions of protons and neutrons are similar for the 
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E [MeV] E [MeV] 

Fig. 3 Fusion cross section for 16 0+ 16 and 16 0+ 18 as a function of the incident energy 
in the center-of-mass frame obtained by adiabatic and sudden potentials for (a) the MV1' and 
(b) D1S interactions. Solid and dotted lines are obtained by sudden and adiabatic 16 0+ 16 
potentials, respectively, and short-dashed and dot-dashed lines are obtained by sudden and 
adiabatic 16 0+ 18 potentials, respectively. The circles and squares denote experimental 
values for 16 0+ 16 and 16 0+ 18 fusion cross sections, respectively, which are taken from 
Ref.[13]. 

16 core, is also shown in the figure. The valence neutron density in 18 has three peaks, 
which is a result similar to that of the intrinsic wave functions of the 18 O ground state 
although the 18 ground state has a spherical density distribution in the laboratory frame. 

We calculated the 16 0+ 16,18 fusion cross sections with the adiabatic and sudden poten- 
tials using the potential model code potfus3[5]. Figure 3 shows the 16 0+ 16 ' 18 fusion cross 
sections as functions of incident energy in the center-of-mass frame. The energy dependence 
of the 16 0+ 16,18 sub-barrier cross sections for £ c . m . < 9 MeV, obtained with sudden and 
adiabatic potentials calculated using the MVl' and D1S interactions, has qualitatively simi- 
lar slopes, except for the 16 0+ 16 fusion cross section calculated using the D1S interaction, 
which decreases rapidly for E cm ^ < 8 MeV. The 16 0+ 16 sub-barrier fusion cross sections 
<7 ad( 16 0) an d o"sud( 16 0) obtained with adiabatic and sudden potentials, respectively, are 
similar to each other. The 16 0+ 18 sub-barrier fusion cross section <r a d( 18 0) obtained by 
adiabatic potentials is a few times larger than the 16 0+ 16 sub-barrier fusion cross sec- 
tions cj a d( 16 0) and <t SU( j( 16 0) because of the lower internuclear potentials. However, the 
cross section cr su( j( 18 0), obtained by sudden potentials, is smaller than cr a d( 18 0) but sim- 
ilar to a a d( 16 0) and <7 su d( 16 0). For the MVl' interaction, the cross sections cr a d( 16 0) and 
0"ad( 18 O) agree with the experimental data qualitatively, whereas for the D1S interactions, 
the estimated sub-barrier fusion cross sections are less than the experimental data. 

4. Discussions 

In this section, we discuss the contribution of valence neutrons in enhancing the sub-barrier 
fusion cross sections by analyzing the results obtained for the MVl' interaction, which 
accounts for the measured 16 0+ 16 ' 18 sub-barrier fusion cross sections qualitatively. As 
mentioned previously, the experimental enhancement of 16 0+ 18 sub-barrier fusion cross 
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Fig. 4 Maximum (dotted) and minimum (solid) 16 0+ 18 internuclear potentials for the 
orientation of 18 0, in which 16 and 18 are frozen, as functions of internuclear distance 
for the MVl' interaction. The sudden (short dashed) and adiabatic (dot-dashed) 16 0+ 18 
potentials are also shown. The long-dashed line represents point Coulomb potential. 

sections compared with 16 0+ 16 fusion cross sections is reproduced by the adiabatic poten- 
tials having lower internuclear potential than the sudden potentials. As shown in Fig. 3, the 
16 0+ 16 ' 18 sub-barrier fusion cross sections obtained with adiabatic and sudden potentials 
are related as follows: 



The relation <r su d( 16 0) ~ a a d( 16 0) indicates that the effects of 16 O excitation on sub-barrier 
fusion cross section are negligible. However, for 16 0+ 18 0, sub-barrier fusion cross sections 
obtained with adiabatic potentials, which is significantly larger than those obtained with 
sudden potentials, indicate that the excitation of 18 enhances the 16 0+ 18 fusion cross 
sections. The relation o" suc j( 16 0) ~ c a d( 16 0) ~ c su d( 18 0) implies that when the 18 nucleus 
is frozen, the theory cannot account for the enhancement of the 16 0+ 18 sub-barrier fusion 
cross sections. Since excitation of the 16 O core is minor, the excitation of 18 O is primarily due 
to the excitation of two valence neutrons around the 16 core. In other words, the excitation 
of the two valence neutrons in 18 O enhances the sub-barrier fusion cross sections, whereas 
increasing the number of neutrons without excitation does not have significant contribution. 

To investigate the details of the excitation effects of valence neutrons in 18 O, we dis- 
cuss the alignment and dipole polarization. Since the 18 ground state exhibits an intrinsic 
deformation, one of the possible excitation effects is the alignment of deformed 18 O. Dipole 
polarization can be another excitation effect. The neutrons in 18 may distribute inward 
(toward 16 O) relative to protons because of the isospin dependence of nuclear interactions 
and the Coulomb force, which results in the isovector dipole polarization of 18 0. Analysis 
of the alignment effects and dipole polarization reveals that the alignment effect contributes 
significantly to the enhancement of the 16 0+ 18 sub-barrier fusion cross sections. 

For calculating 16 0+ 18 sudden potentials, we assumed that nuclei 16 and 18 are frozen, 
and the orientation f2 of 18 is averaged. To study the effect of the alignment of the deformed 
18 on 16 0+ 18 potentials, we analyzed them before averaging the 18 orientation Q. At 



^sud( 16 0) ~ a ad ( 16 0) ~ a sud ( 18 0) < a ad ( 18 0). 



(24) 
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Fig. 5 Dipole moments of 16 (dashed) and 18 (solid) in 16 0+ 18 structures obtained 
by varying the energy for the MVl' interaction. 



every internuclear distance, the orientation Q is optimized to obtain the minimum (max- 
imum) 16 0+ 18 energy Vi' 6 q_ 18 q(R,Q), from which the minimum (maximum) 16 0+ 18 
potential is calculated with frozen 16 and 18 0. The results are shown in Fig. 4 as functions 
of the internuclear distance. The 16 0+ 18 sudden potential (after averaging Q) and the adi- 
abatic potential are also shown for comparison. In the R > 6 fm region, the sudden potential 
is similar to the maximum internuclear potential. In the R > 7 fm region, the 16 0+ 18 adi- 
abatic potential is similar to the minimum internuclear potential, in which 16 and 18 
are frozen. In particular, the minimum values of 16 0+ 18 potential give almost the same 
barrier height and thickness as those of the 16 0+ 18 adiabatic potential, at least for fusion 
at i^c.m. > 7 MeV. These results indicate that the alignment of 18 majorly describes the 
difference between the sudden and adiabatic potentials. We conclude that the enhancement 
of the 16 0+ 18 sub-barrier fusion cross section is due to the alignment effect. 

Dipole polarization can be another effect of excitation in 18 O. Here the isovector dipole 
moments M(£l) are defined as, 

M(El) = ^-^(R p -R n ), (25) 

where R p and R ra denote the positions of the centers of mass of protons and neutrons, 
respectively, and are defined in Eqs. (9) and (10). In the 16 0+ 18 system, 18 valence 
neutrons are somewhat polarized because the existence of 16 O results in finite isovector 
dipole moments for 18 0. Figure 5 shows z-components of the dipole moments M Z (E1) of 
16 and 18 in the 16 0+ 18 wave functions as functions of the internuclear distance R, 
where z-axis is a major axis of the total systems. \M x (El)\ and \M y (El)\ are considerably 
smaller than \M Z (E1)\. The 16 and 18 nuclei are located at z < and z > 0, respectively. 
The 16 0+ 18 wave functions are obtained by the d-constraint AMD method using the 
MVl' interaction [Eq. (11)]. In the large-distance region R > 9 fm, M Z (E1) is positive and 
negative for 16 and 18 0, respectively, indicating that protons distribute outward because of 
the Coulomb force. In this region, the dipole moments are comparable to those in 16 0+ 16 0. 
In the R < 9 fm region, M Z (E1) of 18 gradually increases with a decrease in R. M Z (E1) of 
16 becomes positive and gradually increases. In this region, neutrons in 18 are attracted 
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toward 16 O because of nuclear interactions and push neutrons (or pull protons) in 16 0. 
Although the finite dipole moments are observed in the calculations, the deviation of the 
distance between centers of mass of protons in 16 and 18 from the internuclear distance 
is quite small (less than 0.04 fm in the R > 5 fm region). Therefore, the dipole polarizations 
only have a minor effect on the 16 0+ 16 ' 18 internuclear potentials. 

The above discussions regarding the alignment and the dipole polarization of 18 is based 
on a strong-coupling scenario. In a weak-coupling scenario, these results suggest coupling 
with rotational members such as J 71 = 2 + state in the ground-state band contribute to the 
enhancement of sub-barrier fusion cross sections instead of coupling with J n = 1~ states, 
which is consistent with the CC model study [16]. In heavy and well-deformed systems, 
deformation effects to near- and sub-barrier cross sections are discussed using the orientation- 
average of orientation-dependent cross section with sudden potentials [6]. 

5. Conclusions 

We propose a method to obtain adiabatic internuclear potentials via energy variation with 
the intercluster-distance constraint in the AMD framework. The potentials are applied to 
investigate the sub-barrier cross sections of 16 0+ 16 and 16 0+ 18 through a potential 
model. For the MVl' interaction, the theoretical cross sections agree with the experimen- 
tal data, whereas for the D1S interaction, the theoretical cross sections are less than the 
experimental data. Excitation of valence neutrons in 18 enhances sub-barrier fusion cross 
sections. The alignment of deformed 18 is a dominant excitation effect, while dipole polar- 
ization effects are relatively weak. To understand sub-barrier fusion reactions, the details of 
the structural changes should be considered. The present adiabatic internuclear potentials 
work well to describe sub-barrier nuclear fusions qualitatively. 
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